NACA*ftWFA57A02 


RM A57A02 


i 


j 


NACA 

RESEARCH MEMORANDUM 


AN INVESTIGATION OF WING-BODY JUNCTURE INTERFERENCE 

EFFECTS AT TRANSONIC SPEEDS FOR SEVERAL 

SWEPT- WING AND BODY COMBINATIONS 

By John B. McDevitt and Robert A. Taylor 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

May 16, 1957 

Declassified May 29, 1959 



v 


MCA RM A57A02 


NATIONAL ADVISORY COMMITTEE FOR AEROMUTICS 


RESEARCH MEMORANDUM 


AN INVESTIGATION OF WING-BODY JUNCTURE INTERFERENCE 
EFFECTS AT TRANSONIC SPEEDS FOR SEVERAL 
SWEPT-WING AND BODY COMBINATIONS 
By John B. McDevitt and Robert A. Taylor 


SUMMARY 


This report presents the results of an investigation of wing-body 
juncture interference effects at high subsonic and transonic speeds. 

Force characteristics and pressure measurements for various bodies in 
combination with sweptback wings of aspect ratios 3 and 6 were obtained 
through a Mach number range from 0.8 to 1.2. The bodies used consisted 
of basic Sears-Haack bodies of revolution, the basic bodies indented by 
area- rule concepts at design Mach numbers of 1.0 and 1.2, and the latter 
area-rule bodies reshaped without change in cross-sectional areas (by the 
superposition of slender-body multipoles) so that the flow along the wing- 
body juncture chords at high subsonic speeds more closely resembled that 
for the infinite oblique wing. 

The various area- rule applications resulted in large reductions in 
the transonic wave drag. The reshaping of the M = 1 area- rule bodies 
by application of multipole theory did not result in further drag reduc- 
tions. This result is believed to be due to the fact (as indicated by 
both theoretical and measured pressures) that the M = 1 area- rule appli- 
cations for these particular wing-body combinations result in juncture 
pressure distributions at transonic speeds which are quite similar to the 
pressure distribution of the oblique wing in subcritical flow. However, 
the M = 1.2 supersonic area- rule application was found to give unfavor- 
a Lle juncture-chord pressure distributions. The desired corrections to 
the flow were large and could not be fully achieved by a realistic appli- 
cation of multipole theory. A partial adjustment of pressures was 
attempted which resulted in a slight reduction of wave drag near M = 1. 
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INTRODUCTION 


The possibility of employing sweepback for reducing transonic effects 
followed from the early investigations by Gothert and R. T. Jones (refs. 1 
to 3) of the linearized flow about swept wings (see also the work of 
Neumark , ref. 4). However, the full benefit indicated by simple- sweep 
theory at high subsonic speeds could not be obtained in practice, princi- 
pally because the streamwise pressure distribution near the root of the 
swept wing differs from that for the infinite oblique wing in subsonic 
flow. (The full benefit of sweep is achieved if the wing isobars are 
everywhere parallel to the local sweep of the wing . ) 

The unfavorable pressure distribution near the root of the sweptback 
wing can be corrected either by altering the wing geometry near the root 
or by contouring the body in the vicinity of the wing-body junction. The 
possibility of obtaining favorable junction pressure distributions by 
shaping the sides of the body to conform to the general shape of a stream- 
line on the yawed, infinite wing has been reported on in references 5 
to 7. The design of wing-body junctions by the use of a quasi- cylindrical 
theory involving ring vortices in linearized subsonic flow was developed 
by Kuchemann , references 8 and 9 . (in the Kuchemann method the pressure 
distribution at the wing-body juncture is assumed to be equal to that for 
the center line of the sweptback wing.) A more exact theory for deter- 
mining the linearized, subsonic flow about wing-body combinations has been 
developed recently (ref. 10). 

Although most of the early investigations of the swept-wing problem 
(refs. 5 to 9 and 11 to l4) suffered from the limited Mach number range 
made available by the conventional wind tunnel, significant reductions in 
drag were obtained at high subsonic speeds by the use of body contouring. 
Using rocket-powered models. Pepper (ref. 7) was able to conduct an inves- 
tigation throughout the transonic speed range and found that large drag 
reductions could be obtained by indenting the body so as to approximate 
the streamline flow over the yawed, infinite wing. 

Soon after the completion of the first transonic wind tunnel at the 
Langley Aeronautical Laboratory, the transonic area-rule concept for 
slender- wing-body combinations was discovered and verified by the experi- 
ments of Whitcomb (ref. 15). The extension of this concept to include 
supersonic design Mach numbers was made by R. T. Jones, reference l6 (see 
also ref. 17). Following the discovery of the area rule a re-examination 
of the swept-wing problem appeared desirable. An experimental investiga- 
tion (ref. 18) indicated that both the area- rule concept and Kuchemann' s 
design method result in large drag reductions at transonic speeds for 
swept-wing and body combinations . Furthermore , it was found that the 
area- rule indentation for this particular swept wing resulted in junction 
pressures similar to those for the oblique wing in subsonic flow. Actually 
the body shapes were quite similar, the major difference being in the body 




NACA RM A57A02 


3 


contouring downstream of the junction chord. Thus it appears that for 
swept wings , somewhat similar body contours result from applications of 
the area rule and from the method of KUchemann. 

The possibility of combining KUchemann’s concept of adjusting juncture 
pressure distributions and area-rule concepts motivated the investigation 
of the present report. The two concepts can be considered simultaneously 
by first prescribing the axial distribution of cross-sectional areas 
according to the area rule and then using slender-body multipoles (which 
alter the body shape without change in cross-sectional area) to adjust 
the pressure distribution at the wing-body juncture to conform to those 
for an infinite yawed wing. Several investigations somewhat similar to 
the investigation of the present report have been conducted recently. In 
references 19 to 21 the possibility of shaping the body of sweptback wing 
and body combinations so as to combine the curvature of the streamline 
over the oblique wing with the longitudinal area distribution obtained by 
area-rule concepts has been investigated for both nonlifting and lifting 
cases . 


NOTATION 

The primary symbols used in this report are defined as follows: 

A wing aspect ratio 

A 0 (x) strength of axial sources 

strength of axial multipole distributions according to slender- 
b n (x) J body theory 

b wing span 

Cp drag coefficient, 

Ct lift coefficient, 

L * qS 

pitching-moment coefficient. Pitching moment about c/4 

qSc 

C p pressure coefficient, local pressure minus free-stream static 

pressure divided by q 


c 


local wing chord 
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M.A.C 


} 


mean aerodynamic chord, 


/ 


b/2 


c 2 dy 


■>b/2 


c dy 


l 

M 

Me, 


M. 


des 


reference chord near wing-body junction (chord through the point 
of intersection of basic or area-rule bodies and wing leading 
edge) 

chord at wing center line 
chord at wing tip 

body length (distance from nose to theoretical point of closure) 
free-stream Mach number 
critical Mach number 
design Mach number 

1 2 

free-stream dynamic pressure, -pU^ 


R 

r 

S 

t 

c 

Uoo 

u,v,w 


V 


(r) 


x,y,z 


i 
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body radius 

polar coordinate in y,z plane 
area of wing plan form 

wing thickness-to-chord ratio 
free-stream velocity 

perturbation velocities normalized by division by the free-stream 
velocity 

radial component of perturbation velocity in yz plane normalized 
by division by the free-stream velocity 

Cartesian coordinates 

Cartesian coordinates normalized by division by the juncture 
chord Cj 

Cartesian coordinate with origin at the leading edge of the 

reference chord c . 
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a angle of attack 

p 7im 2 -i| 

9 polar angle in y,z plane 

C T 

A wing plan form taper ratio, — 

C R 

p free-stream density 

cp perturbation velocity potential 

A angle of sweep, positive when swept back 

( ) variable of integration 

( ) ' first derivative with respect to the free-stream direction 
( ) second derivative with respect to the free-stream direction 

Subscripts 

a reference or starting point 

j condition along wing-body junction 

LE wing leading edge 

n order of multipole 

w wing 

2 axial quadripole 


APPARATUS AND MODEIS 
Apparatus 


The tests were conducted in the Ames l4-foot transonic wind tunnel 
which is equipped with a perforated test section permitting continuous 
operation from subsonic to low supersonic speeds . 
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The models were mounted on a sting support as shown in figures 1 
and 2. The normal and chord forces and the pitching moment were measured 
by a strain-gage balance enclosed within the model. Multiple-tube mercury 
manometers, connected to pressure orifices in the model by flexible tubing, 
were photographed to provide records of the pressure distributions on the 


Plan-form details of the various models are presented in figure 3- 
For all the wing-body combinations the center lines of the bodies were 
located in the chord plane of the wing and all bodies were truncated to 
permit mounting on the sting. The body contour details in the vicinity 
of the wing are presented in figures 4 to 6. 

The wing used for model series I (fig* 3( a )) h a< i an aspect ratio 
of 6 a taper ratio of 0.4, and NACA 64A008 sections in the streamwise 
direction with the 40- percent- chord line swept back 40 . The wing used 
for model series II and III had an aspect ratio of 3, a taper ratio of 0.- 
and NACA 64A006 sections perpendicular to the 25 -percent- chord line. The 
wing leading edge was swept back 45°. 

Models I-A and II-A.- The bodies used for models I-Aand II-A were 
bodies of revolution shaped in accordance with the Sears-Haack formula 


The fineness ratio (based on the theoretical length of the body to closure 
and the maximum body diameter) was 11 for body I-A and 12.5 for body II-A. 

Models I -B, II-B, and III-B .- The bodies of revolution used for models 
I-B and II-B were indented according to the M = 1.0 area rule such that 
the wing-body combinations had axial distributions of cross-sectional areas 
equivalent to the basic Sears-Haack bodies of models I-A and II-A, respec- 
tively. (A convenient method for calculating the indented body radii is 
presented in Appendix A of this report.) 

The body of revolution used for model III-B was indented according 
to the supersonic slender-body area mile for a design Mach number of 1.2 
with the area distribution of the Sears-Haack body of model II-A used as 
the desired optimum. (See ref. l6 or 22 for a discussion of the supersonic 
slender-body area rule . ) 


models 


Models 
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Models I-C, II-C, and III-C .- The body shapes for models I-C, II-C, 
and III-C were obtained by altering, without change in cross-sectional 
areas, the shapes of bodies I-B, II-B, and III-B, respectively, by appli- 
cations of s lender-body quadripole theory so as to adjust the flow along 
the junction chords to agree more closely with that for the oblique wings 
in subsonic flow (design Mach number of 0-9). The theoretical method used 
for obtaining the body shapes of models I-C, II-C, and III-C is described 
in detail in Appendix B of this report. 

The calculated pressure distributions (using the wing-body theory 
of ref. 10) along the junction chords of models I-B, II-B, and III-B are 
compared with those of the oblique wing in part (a) of figures 7 to 9 . 

The desired adjustment of juncture pressures and the actual corrections 
considered in the design of bodies I-C, II-C, and III-C are presented 
in part (b) of figures 7 to 9* The computed radii for model II-C near 
9 = 90 ° and 270 ° (top and bottom of fuselage) for body stations 5 h to 68 
(see fig. 5 (h)) were arbitrarily increased a moderate amount to the values 
shown in figure 5 (b), 9 = 90 °, so as to permit installation of the balance 
mechanism. Because of the large differences in pressures indicated in fig- 
ure 9 (a) , the full correction of the junction flow by the use of quadripole 
theory was not possible for model III-C. However, a partial correction was 
arbitrarily chosen as shown in figure 9 . 

In the design of these bodies it was found that the starting, or 
reference, point for the front and rear superposition of the axial quad- 
ripoles must be chosen with care in order to keep the resulting body dis- 
tortions of a realistic magnitude. The various starting points are indi- 
cated in part (b) of figures 7 to 9- Arbitrary fairings consistent with 
slender-body concepts were used to terminate the quadripole distortions 
of the body shapes fore and aft of the junction chord. 


TESTS AM) PROCEDURE 


The models were tested through a Mach number range from 0.8 to 1.2. 

The Reynolds number (based on the wing mean aerodynamic chord and the aver- 
age temperature of the air) was approximately 4xl0 6 for model series I and 
approximately 7xl0 6 for model series II and III (fig. 10). 

Tunnel-boundary- interference corrections were not applied to the 
data. These effects at subsonic speeds are considered to be negligible 
due to the perforated test section. For models of the size employed in 
the present investigation, the influence of the reflected waves on model 
characteristics is known to be small and confined to the Mach number range 
from 1.00 to about 1 . 15 . 

The drag data have been corrected for an interaction in the balance 
mechanism of normal force on chord force and have been adjusted to 
represent free-stream static pressure at the model base. 
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FORCE STUDY RESULTS 


The aerodynamic characteristics of the various wing-body combinations 
of model series I (swept wing of aspect ratio 6) and model series II and 
III (swept wing of aspect ratio 3) are, for convenience, presented sepa- 
rately. Some results of body- alone tests conclude the presentation of 
force data. 


Model Series I 


The drag, lift, and pitching-moment data are shown in figure 11 for 
the three configurations utilizing the swept wing of aspect ratio 6. (See 
figs. 3(a) and h for design details of the various models.) Cross plots, 
which summarize the more important force characteristics, are presented in 
figures 12 to l4. 

Drag.- The variation of drag coefficient with Mach number at zero 
lift and at C L = 0.2 is shown in figure 12. The transonic drag charac- 
teristics of the two models having smooth axial distributions of cross- 
sectional areas in accordance with the M = 1 area rule (models I-B and 
I-C) demonstrate the fundamental importance of the area- rule concept. 

Wear M = 1 the area-rule applications resulted in reductions of wave 
drag, when compared with the basic model I-A, of approximately 70 percent. 
With increasing supersonic Mach numbers, however, the favorable interfer- 
ence effects of the area- rule indentations decrease and eventually 
disappear near M = 1.2. 

The application of multipole theory in the design of model I-C (which 
altered the body cross-sectional shapes without change in the axial distri- 
bution of cross-sectional areas required by the area rule) to improve the 
flow near the wing-body juncture at high subsonic speeds did not appreci- 
ably affect the wave drag at transonic speeds. This result is believed to 
be due to the fact that the M = 1 area- rule application in this particular 
case results in wing-body- juncture pressures at transonic speeds which are 
quite similar to the pressure distribution of the oblique wing in subsonic 
flow. 


The variation of maximum lift-drag ratio with Mach number is shown 
in figure 13- Considerably higher values of maximum lift- drag ratios were 
obtained at high subsonic and transonic speeds as a result of the 
application of the area- rule concept. 

Lift.- The variation of lift- curve slope (evaluated at zero lift) 
with Mach number is shown in figure l4. The two area- rule configurations 
(models I-B and I-C) had, within the accuracy of the experimental data, 
essentially the same values of lift-curve slope throughout the Mach number 
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range. In general, the application of the area- rule concept resulted in 
higher values of lift- curve slope and in a somewhat smoother variation of 
lift-curve slope with Mach number at transonic speeds. This result is a 
direct consequence of the reduced shock losses that invariably accompany 
area- rule modifications for improving the axial distribution of cross- 
sectional areas . 


Model Series II and III 


The drag, lift, and pitching-moment data are shown in figures 15 
and 16 for the five configurations employing the swept wing of aspect 
ratio 3- (See figs. 3(b), 5, and 6 for design details.) Cross plots, 
which summarize the more important force characteristics of these models, 
are presented in figures 17 to 19. 

Drag.- The variation with Mach number of drag coefficient at zero 
lift and at C^ = 0.2 is shown in figure 17(a) for the basic configura- 
tion, model II-A, and for the two configurations, models II-B and II-C, 
having cross-sectional area distributions in accordance with the M = 1 
area rule . The use of the area rule resulted in large reductions of drag 
throughout the transonic speed range. 

The application of a suitable axial distribution of multipoles in 
conjunction with the M = 1 area mile so as to improve the flow along 
the wing-body juncture did not appreciably affect the drag of model II-C. 
As in the previous case for model I-C, this result is believed to be due 
to the fact that the application of the M = 1 area rule for this partic- 
ular wing-body combination results in favorable juncture pressures which 
are difficult to improve . 

In figure 17(b) measured drag characteristics of model II-A are 
repeated with the results for models III-B and III-C, which had area dis- 
tributions according to the supersonic slender-body area rule for a design 
Mach number of M = 1.2. When compared with the previous drag results for 
the M = 1 area- rule designs of model series II, the M = 1.2 area rule is 
found to provide the expected greater wave-drag reductions near the design 
Mach number but is, approximately, only one half as effective near sonic 
Mach number. Theoretical calculations at M = 0. 9, as well as measured 
juncture pressures, indicate that the M = 1.2 area-rule model (model 
III-B) has undesirable flow conditions along the junction chord at tran- 
sonic speeds . The desired correction of the juncture flow was too large 
to be achieved by a realistic application of multipole theory. However, 
a partial adjustment of juncture pressures was attempted by the design of 
model III-C which resulted in a slight reduction of drag near M = 1. The 
variations of maximum lift-drag ratios with Mach number are shown in 
figure 18. 
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Lift.- The variations of lift-curve slope (evaluated at zero lift) 
with Mach number are shown in figure 19 . In general, the various area- 
rule modifications resulted in higher values of lift-curve slope through- 
out the transonic speed range. By inspection of the drag and lift varia- 
tions shown in figures 17 and 19, it is apparent that increases in lift- 
curve slope are directly connected with the relative effectiveness of the 
various area-rule modifications in reducing the transonic wave drag. 


Body-Alone Tests 


The variations of zero-lift drag coefficient with Mach number are 
presented in figure 20 for the various bodies used in this experimental 
investigation. The drag coefficients for body series I and II are based 
on the respective wing areas of model series I and II. 

Bodies I-A and II-A are Sears-Haack bodies of revolution having 
fineness ratios of 11 and 12.5, respectively. Bodies I-B and I-C have 
identical cross-sectional areas with body I-C distorted by an application 
of quadripoles. Bodies II-B and II-C have identical cross-sectional areas 
with body II-C distorted by multipoles; bodies III-B and III-C can be 
described similarly. With the above description of the bodies in mind it 
is apparent that the drag variations shown in figure 20 indicate that the 
wave drag depends primarily on the axial distribution of cross-sectional 
areas. It can also be shown by theoretical considerations, using slender- 
body concepts, that the superposition of axial multipoles does not con- 
tribute to the wave drag (ref. 23), provided that the distribution function 
for the multipole strength, a n (x), terminates with zero first and second 
derivatives . 


PRESSURE STUDY RESULTS 


Surface pressure distributions for the various wing-body combinations 
were measured at zero lift. The experimental results are, for convenience, 
presented separately for model series I (swept wing of aspect ratio 6) and 
model series II and II (swept wing of aspect ratio 3). 


Model Series I 


The measured flow over the wing surface at various Mach numbers is 
illustrated by the isobar patterns presented in figure 21. The locations 
of the isobars were determined from the faired chordwise pressure data 
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obtained at nine lateral stations along the wing semispan. The measured 
pressure distributions along the body surfaces are presented in figure 22 
for two body azmuthal locations, 1 0=0° and 0 = 90°. 

Wing isobars .- The isobars shown in figure 21 demonstrate that the 
pressure field of the wing is decisively influenced at near-sonic speeds 
by the adjacent body shape. The unfavorable distortion of the wing pres- 
sure field near the basic body of model I- A is evident throughout the 
Mach number range. The M = 1 area- rule modification of model I-B cor- 
rected the unfavorable pressure field near the wing-body juncture in a 
very satisfactory manner (that is, the flow is in close agreement with 
that for the oblique wing) . In particular, it should be noted that the 
streamwise location of the point of minimum pressure (indicated by the 
dashed lines of fig. 21) was shifted to a forward location at the root 
chord. The quadripole distortion used in the design of model I-C did 
not significantly alter the isobar patterns. 

Body pressure distributions .- The importance of the body shape is 
again evident from the pressure distributions presented in figure 22. It 
is apparent that the favorable pressure distributions resulting from the 
area- rule modification of model I-B would be difficult to improve in order 
to achieve a further reduction in wave drag near M = 1. 

Comparison of calculated and measured juncture pressures . - The 
theoretical and experimental pressure distributions near the wing-body 
junctions of model series I at zero lift are shown in figure 23 for Mach 
numbers of O.85 and 0.90. The theoretical results were obtained by the 
use of the linearized flow theory for wing-body combinations described 
in reference 10. 

The measured pressures at the innermost row of wing orifices 
(y = 0.139 b/2) are also included in figure 23 since it was felt that a 
comparison of theory and experiment at the wing-body junction would suffer 
somewhat from boundary- layer effects in the corners of the wing-body 
juncture. In general, the theory and experiment are in reasonable 
agreement . 


Model Series II and III 


The measured pressure distributions at various Mach numbers for the 
models of series II and III at zero lift are presented in figures 2k to 29 . 

x At the wing-body juncture the pressure orifices along the sides of 
the bodies were located as close as possible to the upper surface of the 
wing. The designation 0=0° for the side row of orifices is used here 
for convenience, although, at the wing-body juncture, the orifice loca- 
tions differed slightly from 0=0° because of the finite thickness of 
the wing. 
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Wing isobars for models II-A, II-B, and II-C .- The wing pressure 
distributions for models II-A, II-B, and II-C are indicated by the isobar 
patterns shown in figure 24. The unfavorable pressure field near the 
wing-body juncture of model II-A is evident throughout the Mach number 
range. The M = 1 area- mole modification of model II-B reduced the gen- 
eral level of the disturbance field over the entire wing for Mach numbers 
up to 1.10 and improved the flow near the juncture throughout the Mach 
number range. 

Body pressure distributions for models II-A, II-B, and II-C .- The 
measured body pressure distributions for models II-A, II-B, and II-C are 
presented in figure 25. The M = 1 area- rule modification resulted in 
large changes in the body pressures which would be, as for the previous 
case of model I-B, difficult to improve upon in order to achieve further 
drag reductions near M = 1. 

Comparison of calculated and measured juncture pressures for models 

II- A, II-B, and II-C .- The theoretical and experimental pressure distri- 
butions near the wing- body junctions of models II-A, II-B, and II-C at 
zero lift are presented in figure 26 for Mach numbers of O .85 and 0.90. 

In general, the agreement between theory and experiment is good for 
models II-A and II-B. 

Wing isobars for models III-B and III-C .- Wing isobars for models 

III- B and III-C are presented in figure 27. For the M = 1.2 area- rule- 
indented model III-B a local region of positive pressure coefficients 
occurs near the mid- chord of the wing-body juncture where the indented 
body radius is a minimum (see fig. 6(a)). As stated previously, the body 
used in model III-C was the result of an attempt to correct partially the 
junction flow at M = 0.90 by an application of quadripole theory. 

Body pressure distributions for models III-B and III-C .- The measured 
body pressure distributions for models III-B and III-C at zero lift are 
presented in figure 28. Whereas the previous applications of the M = 1 
area rule resulted in very satisfactory pressure distributions at the 
junction chords, it is evident from the pressure distributions for model 
III-B shown in figure 28 that the M = 1.2 area rule is unsatisfactory in 
this respect. The partial application of quadripole theory used in the 
design of model III-C improved considerably the pressure distribution 
along the junction chord. However, the flow along the top portion of the 
body surface was adversely affected by the quadripole distortion. 

Comparison of calculated and measured junction pressures for models 
III-B and III-C.- The theoretical and experimental pressure distributions 
near the wing-body junctions of models III-B and III-C at zero lift are 
presented in figure 29 for Mach numbers of O .85 and 0.90. The agreement 
between theory and experiment is good for model III-B but not so good for 
model III-C where the theoretical pressure coefficient contained large 
contributions from the quadripole distortions of this particular body 
shape . 
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Body-Alone Tests 


The measured surface pressures for the basic bodies (Sears-Haack 
bodies of revolution having fineness ratios of 11 and 12.5) tested alone 
(wings removed) are presented in figures 30 and 31. The pressure distri- 
butions for the two bodies indented by applications of the M = 1 area 
rule are presented in figures 32 and 33. The pressure distributions for 
the body indented by the M = 1.2 slender-body area rule are presented 
in figure 3^* 

It is interesting to note that, for a given indented body shape, the 
streamwise pressure distribution does not change much throughout the Mach 
number range tested. 


CONCLUDING REMARKS 


The theoretical calculations at high subsonic Mach numbers as well 
as the experimental measurements at transonic speeds indicate that the 
M = 1 area- rule designs for sweptback wings result in wing-body- juncture 
pressures quite similar to the streamwise pressure distribution of the 
oblique wing in subcritical flow. As a result it can be expected that, 
for the particular case of the sweptback wing, the area rule will give 
results similar to those obtained by Kuchemann’s method and hence increase 
the critical Mach number as well as provide reductions in transonic wave 
drag. 


In the present experimental investigation it was found that the 
quadripole distortions of the body shapes had comparatively little effect 
on the transonic wave drag when compared with the large effects resulting 
from changes in body cross-sectional areas (which is equivalent to speci- 
fying changes in the axial distribution of sources according to linearized 
or slender-body theory) . It also appears from this investigation that 
significant improvements in M = 1 area-rule applications for transonic 
speeds are improbable for swept wings having at least a moderate amount 
of sweep but that further work is needed to improve juncture pressures 
for wing-body combinations designed by applications of supersonic area 
rules. In recent work by Lomax (ref. 2b) it was shown that the supersonic 
extension of Whitcomb’s area rule is an approximation for the correct rule 
of linearized theory which requires knowledge of model surface pressure 
distributions . Some indication of the success of a particular supersonic 
area-rule application at transonic speeds may be obtained by an inspection 
of the juncture pressure distribution at high subsonic speeds . 

It is probable that a more exact understanding of the sonic and 
supersonic flows about wing-body combinations must be achieved before 
methods for obtaining optimum body shapes can be f ully developed. 
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Attention should also be directed toward the more general problem of 
determining wing geometry such that the wave-drag reductions to be 
obtained shall be as complete as possible. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Jan. 2, 1957 
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APPENDIX A 

BODY- INDENTATION FORMULAS FOR SWEPT- WING AND 
BODY COMBINATIONS INDENTED ACCORDING 
TO THE M = 1 AREA RULE 


Formulas for obtaining the body shape of configurations indented 
according to the M = 1 area rule are easily obtained in the case of the 
swept wing if the wing has straight leading and trailing edges and tip 
chords parallel to the body axis. The area-rule concept suggests that 
the axial distribution of cross-sectional area for the wing- body combina- 
tion be optimized , that is, made equivalent to that for a smooth body 
having low wave drag. The usual procedure is to remove from the smooth 
body volume the exposed volume of both wing panels such that the cross- 
sectional areas in planes normal to the body axis for the indented con- 
figuration are equal to the cross-sectional areas of the smooth body. A 
convenient method for determining the indented radii follows . 


WINGS OF CONSTANT THICKNESS TO CHORD RATIO 


Consider the wing-body 
combination shown in sketch (a) 
where, for convenience, all dimen- 
sions have been normalized by divi- 
sion by the junction chord c*. If 

J 

the integral representation for the 
wing cross-sectional area in the 
plane normal to the body axis is 
replaced by an equivalent integra- 
tion performed in the streamwise 
direction, then the indented body 
radii, Rj, can be obtained by use 
of the following equation: 


Note All distances dimensionless with respect to c 



Body axis 


Sketch (a) 



V R 

cot A 


cot A 1 c 1 s 2Z 1 d§ 1 


(Al) 
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where 2Z X is the wing thickness measured along the junction chord and the 
local chord c x is given by 


c ± = l-K(|j-| 1 )cot A x 


Bci 

<3t) 


1-C n 


(b/2)-R £ 




(A2) 


For all practical cases, equation (Al) can be approximated by the 
following relationship 

\ (R 2 -Rj 2 ) - P(fi J+-TT" F ’(t j) (A3) 

2 ± J cot Aj J 

where 

FUj) = J j cot A 1 c 1 2 2Z 1 d^ 1 

o 


= f .^K^.fp+K 2 ^ . 2 f 3 +2Kf 4 -2K 2 ^.f 5 +K e f 6 {kb) 

J ^ J • J 


F’(|.) «f x '(6j) 


(A5) 


and the f functions and corresponding first derivatives are given by 




- 2Z x cot A 1 
fg'da.) = 2Z x cot 2 A x 
f 3 '(li) = 2Z 1 cot 3 A 1 
f 4 '(5i) = 2Z x | x cot 2 A x 
f 5 *(6 x ) = 2Z 1 5 1 cot 3 A 1 
fe'Ui) = 2Z 1 | 1 2 cot 3 A 1 




(a6) 
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(A7) 


The evaluation of the integrals of equation (A7) usually requires 
numerical or graphical techniques. A convenient formula for obtaining 
the variation of A x with £ is 


tan A x = tan A-^-K^ 


(A8) 


The f i and fb functions for the wing used in model series II of this 
report are shown in figure 35 . 

Solving equation (A3) for Rj gives the desired body- indentation 
f ormula 


RjU.) - 




rt cot A 


j 


.rt cot Aj_ 


+R 2 (!<) 


2F(I<) 2R a fi*(S 1 h X/2 


it cot A 




0 = h = e T 


(A9) 


The tip effect can be accounted for be defining the f functions to 
be 


f i ( f i)g I j'if ’ ^T = = 


A - • 15 !j SfT 

J 


f i _ 5j-5 T ’ 1 < 


(Aio) 


(All) 

(A12) 


18 


NACA RM A57A02 


and noting that in all cases 

f i > = 0 , £j £ 1 ( A1 3) 


WINGS OF VARYING THICKNESS TO CHORD RATIO 


If the wing thickness to chord ratio varies linearly from root to 
tip, the preceding method for finding the indented body radii can he used 
by first replacing Zj_ in equation (Al) by the following 


. hd c j Zl 

a-\. i-?-c V 


where A c is the ratio of tip chord to junction chord and A t is the 
ratio of tip thickness to chord ratio to junction thickness to chord 
ratio. The corresponding f' and f functions are then evaluated and 
used in conjunction with equations (A9) to (A13). 
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APPENDIX B 

THE CALCULATION OF THE BODY SHAPES 
FOR MODE IS I-C, II-C, AND III-C 


The plan form and. body- contouring details of models I-C, II-C, and 
III-C are presented in figures 3> Mb), 5(b), and 6(b). The bodies used 
for these models were obtained by altering, without change in cross- 
sectional areas, the body shapes of models I-B, II-B, and III-B, respec- 
tively, by the application of axial multipoles so as to adjust the flow 
along the junction chord to agree identically with the flow over similar 
oblique wing sections for a design Mach number of 0.90. 


A coordinate system made dimensionless with respect to the junction 
chord will be used, see sketch (b) . The notation of reference 10 will be 
employed and the reader is referred to reference 10 for further details 
concerning the use of planar sources and 
axial multipoles to describe the sub- 
sonic flow about wing-body combinations. 


The first step toward a solution 
of the inverse problem described here 
is the calculation of the juncture 
pressure distribution due to the wing- 
alone sources and body axial sources 


(0 = 0°, M = 0.90), 


= -2(u w+ u 0 )- 


n=o 


-.2 


R f +Av 


(r) 


wJ 


where u. 


L w 


(Bl) 

is the streamwise perturba- 
tion velocity induced by the wing-alone 
sources, u Q is the streamwise pertur- 
bation velocity induced by the axial 
sources, R f is the slope of the area- 
rule body, and the incremental radial 
velocity induced by the wing is defined 
as 



Sketch (b) 


Av 


(r) 


= V 


w 


(r) 


w 


V 


(r)J 

WJ av 


where 


r (r) 


is the radial velocity induced by the wing and 


w 


) W J 

W— 1 c 


is the average radial- induced velocity (see ref. 10). 
the axial sources to be used in the calculation of u<- 


V (r). . 

av 

The strengths of 
are given by 


J 


r (r) 


w- 


av 


20 


NACA RM A57A02 


The next step involves the placement along the body axis of a suitable 
axial distribution of multipoles such that the juncture pressure distribu- 
tion will be identical with that for the infinite oblique wing. The incre- 
mental adjustment of juncture pressures to be accomplished by the axial 
multipoles is therefore 


AC p = C p 

j oblique wing 



n=o 


(B2) 


Since the influence on the pressure field decreases rapidly with the 
order of the multipole (ref. 10), it appears reasonable to consider only 
the multipole of order n = 2. The desired adjustment of the pressure 
coefficient is then (0 = 0°) 



- 21 * 2 - 


R' 



(B3) 


where R f is the slope of the desired body radius. 

In the inverse problem described here, the axial quadripoles are used 
to adjust locally the pressure field but are not used to cancel any part 
of the wing-body interference velocities; consequently, the desired shape 
of the body is given by 


where 


R(|j,0) = R*U j ,0)+AR 2 (| j ,0) 


R*(8j,e) 



(B4) 


(B5) 


(The lower limit of the above integral is arbitrarily taken at £ = -0.5 
since the influence of the wing is negligible forward of this point.) 

The parameter AR 2 in equation (B^) denotes the incremental change in 
body radius due to the axial multipole. 

If use is made of the definitions of equations (b 4) and (B5), 
equation (B3) can be written as (0=0°) 



= -2u 2 - (AR 2 1 ) 2 -2AR 2 t R* f 


(b6) 
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where 


ARg ' = R'-R*’ (B7) 

The desired body shape is to be obtained by solution of equation (B6) . 
First, however, it is necessary to relate the streamwise perturbation u 2 
to the parameter AR 2 . This is easily achieved if the following slender- 
body approximation for multipole theory (ref. 25) is used: 


<P n (x,r,e) 


a n (x)cos n0 b n (x)sin n 6 
“+ 


„n 


„n 


(b8) 


The streamwise and radial perturbation velocities for the quadripole are 


u 2 (x,r,0) = 


a 2 ? (x)cos 2 9 


(B9) 


v (x,r,0) = 


-2a 2 (x)cos 20 


(BIO) 


The axial strength distribution for the quadripole can be related to 
body geometry at 0=0° as follows: 


-2a„(x) 


(r). 


= R'-R*' 


R' 


= AR 2 * 


a 2 (x) 


-R 3 AR 2 ' 
2 


(Bll) 


a 2 * (x) 
R 2 


-R£R 2 " 


-3R *AR 2 * 


(B12) 


2 


2 
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It is now possible to express equation (b 6) as 

AC P = R*£R 2 " +AR 2 AR 2 " +R* 'AR ’ +2 (AR 2 ’ f (B13) 

i 


which can be approximated by 


ACpj 

R* 


R* ’ARo 1 

£R 2 " + ~ 

R* 


(BlM 


(Eq. (Bl4) indicates that the pressure adjustment duetto an axial 
distribution of quadripoles for bodies having small R' is directly 
proportional to the second derivative of the radial modification.) 

For the application described here (the design of models I-C, II-C, 
and III-C) equation (Bl4) was solved by an iteration procedure starting 
with the approximation 


ar 2 * 




R* 





AR 2 'd£ 


where 


tj -'5 a 


is the starting (or reference) point and 


(B15) 


(Bl6) 


AR 2 (l a ) = AR 2 ‘(| a ) = 0 

The starting points (see figs. 7 "to 9) were found to be somewhat critical 
if the body shapes were to remain practical. Solutions of equation (Bl4) 
were obtained along the wing-body junctures and arbitrary fairings were 
used to terminate the body modifications fore and aft of the junction 
chords . 

The solution for AR 2 at 9 = 0 is equivalent to finding the strength 
distribution for the quadripoles. At finite 9 values, AR 2 can be 
obtained by an iteration of the following: 
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r (r), 


■2a 2 cos 2 0 

¥ 


(B17) 


R = R*+ 



I 


a 



(B18) 


An alternate iteration procedure for obtaining the desired body shape 
is obtained if equation (B17) is written in the form 


R'-R* ' 


-2a 2 cos 20 
R 3 


or 


R 3 R* = R 3 R* -2a 2 cos 20 


(B19) 


and finally 


R 4 (Sj) ■ 




I. 


R 4 (i a )+^ / ^ R 3 R* ,( 3-i _ 8 cos 29 J J a^g^dl 


1/4 


(B20) 


As a first approximation, R in the right-hand side of equation (B20) 
can be replaced by 



R = R*-2 cos 20 


(B21) 
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Figure I- General arrangement of the test section of the Ames 14 -foot transonic wind tunnel. 
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Figure 2.- A sting -mounted model in the Ames 1^-foot transonic wind tunnel 
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(a) Model series I. 

Figure 3.- Plan-form details of the models. 
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(b) Model series II and HE. 
Figure 3 - Concluded. 
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(a) Models I-A and I-B. 
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Body station, inches 

(b) Model I-C. 

Figure 4.- Body contour details in the vicinity of the wing for model series I. 
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Body station, inches 
(b) Model IT-C. 

Figure 5.- Body contour details in the vicinity of the wing for model series IT. 
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Body radius, inches 






Body station, inches 
(b) Model HT-C. 


Figure 6.- Body contour details in the vicinity of the wing for model series IH. 
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Junction-chord station, £j 

(a) Calculated pressure distributions at the wing-body juncture 
of model I-B and for the oblique wing. 
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-4 -.2 0 .2 .4 .6 .8 1.0 1.2 1.4 

Junction-chord station, £j 

(b) Incremental adjustment of junction pressure required 
by the axial quadripoles of body I-C. 

Figure 7.- Theoretical calculations (M d e« = 0.90) considered in the body 

design of model I-C . 
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-.4 -.2 0 .2 .4 .6 .8 1.0 1.2 1.4 

Junction-chord station, £j 

(a) Calculated pressure distributions at the wing-body juncture 
of model D-B and for the oblique wing. 



-.4 -.2 0 .2 .4 .6 .8 1.0 1.2 1.4 

Junction -chord station, £j 

(b) Incremental adjustment of junction pressure required 


by the axial quadripoles of body H-C. 

Figure 8.- Theoretical calculations (Md»« = 0.90) considered in the body 

design of model n-C. 
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-4 -2 0 .2 .4 .6 .8 1.0 1.2 1.4 

Junction -chord station, £j 

(a) Calculated pressure distributions at the wing-body juncture 
of model H-B and for the oblique wing. 



Junction -chord station, £j 

(b) Incremental adjustment of junction pressure required 
by the axial quadripoles of body m-C. 


Figure 9.- Theoretical calculations (M d e«=0.90) considered in the body 

design of model DI-C. 
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Mach number, M 


Figure 10.- The variation of average test Reynolds number with Mach number. 
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(a) Drag. 


Figure II - Drag, lift, and pitching-moment characteristics of the models of series I. 
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(a) Drag — Concluded. 


Figure II.- Conti nued. 
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(b) Lift. 


Figure 1 1 - Conti nued. 
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(b) Lift — Concluded. 

Figure II.- Continued. 
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(c) Moment. 

Figure II. - Continued. 
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(c) Moment — Concluded. 


Figure II Conclu ded. 
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Mach number, M 

Figure 12- The variation of drag coefficient with Mach number at zero lift 
and C l = 0.2 for the models of series I. 


Maximum lift-drag ratio, (L/D)max 



Figure 13.- The variation of maximum lift-drag ratio with Mach number for the models 

of series I. 
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Figure 14.- The variation of lift-curve slope with Mach number for the models of series I. 
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(a) Drag. 

Figure 15 - Drag, lift, and pitching-moment characteristics of the models of series!. 
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Figure 15.- Continued. 
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(b) Lift. 

Figure 15 - Continued. 
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Figure I5. - Continued. 
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(c) Moment. 

Figure 15- Continued. 
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(a ) Drag. 


Figure 16 Drag, lift, and pitching-moment characteristics of the models of series III . 
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(a) Drag — Concluded. 
Figure 16.- Continued. 
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Figure 16 .- Continued. 
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(b) Lift — Concluded. 
Figure 16 .- Continued. 
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(c ) Moment. 
Figure 16.- Continued. 
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Figure 16“ Concluded. 
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(a) Models D-A, H-B.and H-C. 

Figure 17- The variation of drag coefficient with Mach number at zero lift 
and C L = 0.2 for the models of series H and series HI. 
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( b ) Models II-A, EI-B, and IH-C. 

Figure 17- Concluded. 
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Figure 18- The variation of maximum lift-drag ratio with Mach number 
for the models of series II and series HI. 
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Figure 19.- The variation of lift-curve slope with Mach number for the models 

of series II and series HI. 
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Figure 20.- The variation of zero-lift drag coefficient ( based on the 
respective wing areas) with Mach number for the various bodies. 
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Figure 21 - Experimental wing isobars for model series lot 

zero lift. 
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Figure 21.- Concluded. 
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(a) M = 0.90 

Figure 22.- Experimental body pressure distributions for model 

series I at zero lift. 
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Figure 22.- Continued. 
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(d) M = 1.02 
Figure 22.- Continued. 
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Figure 22.- Continued. 
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(f) M = l.l 0 

Figure 22.- Concluded. 
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(a) M = 0.85 

Figure 23.- Theoretical and experimental pressure distributions near 
the wing-body junctions of model series I at zero lift. 
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Figure 23.- Concluded. 
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(a) M =0.90 

Figure 24.- Experimental wing isobars for 
models H-A.II-B, and D-C at zero lift. 
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Figure 24.- Continued. 
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Figure 25- Experimental body pressure distributions for models H-A, 

n-B, and H-C at zero lift. 
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(b) M = 0.95 

Figure 25.- Continued. 
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(c) M = 1.00 

Figure 25.- Continued. 
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Figure 25.- Continued. 
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Figure 25 .- Continued. 
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Figure 25.- Concluded. 
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(a) M = 0.85 

Figure 26 - Theoretical and experimental pressure distributions near the 
wing-body junctions of models n-A, n-B, and n-C at zero lift. 
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Figure 26 .- Concluded. 
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Figure 27.- Wing isobars for models 
m-B and m-C at zero lift. 
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Figure 27.- Continued. 
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(d) M = 1.15 

Figure 27.- Concluded. 
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Figure 28- Experimental body pressure distributions for models m-B 

and m-C at zero lift. 
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(b) M = 0.95 

Figure 28.- Continued. 
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Figure 28.- Continued. 
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Figure 28.- Concluded. 
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(a) M = 0.85 

Figure 29- Theoretical and experimental pressure distributions near the 
wing-body junctions of models m~B and HI-C at zero lift. 
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(b) M = 0.9 0 
Figure 29.- Concluded. 
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Figure 30.- Experimental pressures for body I-A (Sears-Haack body of 

fineness ratio II) at zero lift. 
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Figure 31 - Experimental pressures for body H-A (Sears-Haack body of 

fineness ratio 12.5) at zero lift. 


Pressure coefficient, 


NACA RM A5TA02 


10 ^ 




Figure 32.- Experimental pressures for body I-B at zero lift. 
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Figure 33 - Experimental pressures for body II-B at zero lift. 
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Figure 34 - Experimental pressures for body HL-B at zero lift. 
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Figure 35.- The f,-' and f; functions (eqs.A6 andA7) for the 
wing of model series H . 


NACA - Langley Field, Va. 


